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Abstract 

We examine electron and nuclear recoil backgrounds from radioactivity in the ZEPLIN-III dark mat- 
ter experiment at Boulby. The rate of low-energy electron recoils in the liquid xenon WIMP target is 
0.75±0.05 events/kg/day/keV, which represents a 20-fold improvement over the rate observed during the 
first science run. Energy and spatial distributions agree with those predicted by component-level Monte Carlo 
simulations propagating the effects of the radiological contamination measured for materials employed in the 
experiment. Neutron elastic scattering is predicted to yield 3.05±0.5 nuclear recoils with energy 5-50 kcV 
per year, which translates to an expectation of 0.4 events in a 1-year dataset in anti-coincidence with the 
veto detector for realistic signal acceptance. Less obvious background sources are discussed, especially in 
the context of future experiments. These include contamination of scintillation pulses with Cherenkov light 
from Compton electrons and from (3 activity internal to photomultipliers, which can increase the size and 
lower the apparent time constant of the scintillation response. Another challenge is posed by multiple-scatter 
7-rays with one or more vertices in regions that yield no ionisation. If the discrimination power achieved 
in the first run can be replicated, ZEPLIN-III should reach a sensitivity of ^lxl0~ 8 pb-year to the scalar 
WIMP-nucleon elastic cross-section, as originally conceived. 
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1. Introduction 

The ZEPLIN-III direct dark matter search has 
been operating at the Boulby underground labora- 
tory (UK) under a rock overburden of 2,800 m wa- 
ter equivalent. The WIMP target is a xenon emis- 
sion detector [l[ which discriminates between elec- 
tron and nuclear recoils by measuring the relative 
amount of prompt scintillation light and ionisation 
charge extracted from particle interactions in the 
liquid xenon (LXe) phase; the ionisation is trans- 
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duced into a second optical signal via electrolumi- 
nescence in the thin vapour layer above the liquid. 
Both the primary (SI) and secondary (S2) scintilla- 
tion pulses are detected by an array of 31 photomul- 
tipliers immersed in the liquid. Key components of 
the experiment are shown in Figure [TJ for details on 
the design and construction of the system we refer 
the reader to Refs. @, 0, i, H . 

After deployment underground in early 2007, an 
83-day run produced competitive limits for WIMP- 
nucleon scattering cross-sections 0, 0, @|- In this 
first phase the background of the experiment was 
dominated by the photomultiplier array by a signif- 
icant factor - both in electron recoils from 7-rays 
and in nuclear recoils from neutrons. A develop- 
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Figure 1: Schematic diagram of the ZEPLIN-III experiment. 
1: WIMP target (LXe shown in blue); 2: photomultiplicr 
array; 3: xenon chamber; 4: vacuum cryostat; 5: Gd- 
loaded polypropylene; 6: veto scintillator modules; 7-bare 
polypropylene; 8-copper flooring; 9-lead shielding. 



ment programme was thus undertaken with indus- 
try in order to develop new phototubes targeting an 
order of magnitude reduction in radioactivity. After 
upgrade of the photomultipliers and the addition of 
a veto system around the main target [1, [B| , the ex- 
periment returned to taking WIMP search data in 
mid 2010. The first aim of this article is to describe 
the radioactivity backgrounds which determine the 
sensitivity of the rare event search ahead of forth- 
coming results from the second run. 

We highlight also that many background stud- 
ies focus solely on simple event topologies quan- 
tifiable by Monte Carlo simulation, but often over- 
look other important factors. Although the require- 
ment of negligible neutron background drives the 
design of most dark matter experiments, in fact 
this threat materialises only rarely. In two-phase 
xenon detectors, multiple scattering 7-rays can be 
challenging when one of the vertices occurs in a re- 
gion which yields no ionisation; in this instance the 
prompt scintillation signals for the multiple vertices 
are time-coincident and difficult to distinguish, and 
the charge-to-light ratio (S2/S1) is low, more typi- 
cal of nuclear recoils. In another example, when es- 
timating photomultiplicr backgrounds, many stud- 
ies consider only high-energy 7-rays as contribut- 



ing significantly to electron recoil background; we 
present evidence here that /3-induced signals gener- 
ated internally (involving, for example, Chcrcnkov 
emission from the phototube envelopes) should also 
be assessed carefully. 

Besides their sensitivity in the scintillation and 
ionisation channels (see, e.g., Ref. Q for an em- 
phatic demonstration of the latter), two-phase 
xenon detectors hold great promise for next gen- 
eration dark matter searches due to the ability to 
self-shield against external backgrounds. The accu- 
rate reconstruction of the position of interactions 
in three dimensions allows a sacrificial LXe vol- 
ume to be defined which shields an inner 'fiducial' 
mass from most sources of radioactivity. ZEPLIN- 
II was the first WIMP experiment to operate using 
this technology fill , soon followed by XENON10 
at Gran Sassofllf and ZEPLIN-III at Boulby @. 
XENON100 [ljalready benefits from self-shielding 
very significantly, and even more so will the up- 
coming LUX350 experiment at SUSEL [Hj]. Two- 
phase xenon experiments have been at the forefront 
of WIMP sensitivity and systems with tonne-scale 
fiducial mass could probe most of the parameter 
space favoured by constrained minimal supersym- 
metry (cMSSM) and by other extensions to the 
standard model. With this work we wish to sum- 
marise a decade of experience in designing, building 
and operating ZEPLIN-III, and thus inform the de- 
sign of next-generation WIMP experiments. 

This article is organised as follows. Ordi- 
nary, single-vertex electron recoil backgrounds are 
treated in Section ^ where ZEPLIN-III data are 
confronted with an overall prediction built up from 
component-level simulations. In Section [3] we dis- 
cuss event topologies that may degrade the discrim- 
ination power of these instruments, such as 7-ray 
multiple scatters and (3~ background internal to the 
photomultipliers. Having validated U/Th contami- 
nation levels in key components, we present calcula- 
tions of neutron background in Section [4] Finally, 
we summarise our findings in Section [5l with the 
design of future systems in mind. 

2. Single-vertex electron recoil backgrounds 

Most low-energy electron recoils in the LXe re- 
sult from interactions of 7-rays and /3-particles from 
natural radioactivity. We begin by calculating ab- 
solute spectra of deposited energy from dominant 
sources, and compare these with ZEPLIN-III data 
acquired at Boulby in fully-shielded configuration. 
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These signals can be discriminated from nuclear re- 
coils to a large degree; it is reasonable to expect 
that a 'leakage' (mis-identification probability) of 
order 1:10,000 may be achieved for single-site inter- 
actions such as those from 85 Kr /3-particles within 
the LXe (the value eventually achieved in the first 
run for all interactions was 1:7,800 [HI). However, 
more complicated topologies must be considered, 
and these may not be rejected so efficiently. We 
analyse two such types in Section [3] 

The 7-ray background is dominated by primor- 
dial U, Th and 40 K. We assume natural uranium 
(0.956 238 U, 0.044 235 U) and that all chains are in 
secular equilibrium, except where indicated other- 
wise. To determine contamination levels, a large 
number of materials were radio-assayed over the 
lifetime of the project. All HPGc 7-ray spec- 
troscopy was conducted at Boulby with an OR- 
TEC GEM detector with a 2 kg crystal, reach- 
ing down to ~0.5 ppb sensitivity for U and Th. 
Energy-dependent counting efficiencies were calcu- 
lated from Monte Carlo simulations for each test ge- 
ometry. These measurements were complemented 
by commercial mass spectrometry (ICP-MS) for 
several materials. We have also benefited from the 
extensive UKDMC database [H] . 

The 20 most intense 7-rays in each chain and the 
1,461 keV line from 40 K are propagated using the 
GEANT4 Monte Carlo toolkit [l6| to interactions 
in the LXe target. For internal backgrounds the 
ZEPLIN-III simulation model was used [2| ; this was 
complemented by importing a detailed CAD solid 
model for calculations involving the veto, shielding 
and the laboratory rock. 

2.1. ZEPLIN-III measurements 

A 24-day long dataset was analysed to study the 
spatial and energy distribution of electron recoil 
background. This followed a neutron calibration 
by only a few weeks, and so isomeric transition 
(IT) 7-rays from activation of the xenon are still 
visible (mainly from 129m Xe, T!/ 2 =8.88 days, and 
131m Xe, Ti/ 2 =11.8 days). The data were collected 
as part of the WIMP search dataset and, for this 
analysis, the nuclear recoil signal region remained 
blinded up to 40 kcVcc0 However, this has little im- 
pact on the electron recoil population. Events with 



1 LXe energy deposits are expressed in electron-equivalent 
'keVee' units when calibrated by 57 Co 122 keV 7-rays, or in 
'keVnr' when referring to nuclear recoil energy. 
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Figure 2: Electron recoil background in ZEPLIN-III. Data 
are shown in yellow and simulated components are summed 
by the solid blue line, with individual contributions ordered 
by importance at low energy. Top: energy spectrum in fidu- 
cial volume, showing good agreement below ~150 keVee, 
the dynamic range of the dataset (the thin black histogram 
retains only unsaturated events); U/Th and Xe activation 
7-rays are indicated. Middle: radial distribution below 
100 keVee, confirming very good agreement beyond the cen- 
tral channel. Bottom: depth distribution for fiducial region 
(bright yellow) and full depth (hatched yellow); contami- 
nation of the cathode grid (37 mm) with radon progeny left 
over from the first run is clearly visible (but easily fiducialiscd 
away) . 
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single SI and S2 pulses are selected, which allows 
three-dimensional reconstruction of the interaction 
point. A 6.5 kg fiducial region (150 mm radius, 
32 mm deep) is retained for analysis, with no signifi- 
cant quality cuts applied. The energy variable used 
is a linear combination of SI and S2 energy esti- 
mators which exploits the anti-correlation between 
these two response channels (ionisation electrons 
may contribute either to S2, if extracted from the 
track, or to SI, via recombination luminescence). 
This allows for a better measurement of energy than 
cither channel independently. The dynamic range 
of the dataset is limited, owing mainly to the 8- 
bit resolution of the digitisers. Although a dual- 
range acquisition system is used, the large S2 gain 
means that signal truncation becomes significant 
above ~150 keVee in the WIMP-search dataset. 

The measured energy and spatial distributions 
are shown in Figure [5] - along with a component- 
level breakdown from simulated data calculated 
as described below. These contributions are also 
listed in Table I. The measured low-energy back- 
ground is 0.75±0.05 events/kg/day/keVee (here- 
after 'dru', for 'differential rate unit'). Simulated 
and measured spectra are in good agreement up to 
~150 kcV and even beyond, despite the xenon acti- 
vation which is still prominent in this dataset. This 
represents a ^20-fold improvement over the first 
run. The simulation sum predicts 0.85±0.05 dru 
at low energy; as discussed further ahead, this is 
based on measured contamination only, i.e. exclud- 
ing items for which contamination levels could not 
be detected. Aside from this selection, there is no 
scaling of the simulation results in that figure. 

The spatial distributions are also in good agree- 
ment, with two exceptions. The radial trend 
matches, within errors, over most of the fiducial vol- 
ume, but the rate is slightly higher than predicted 
within the radius corresponding to the central pho- 
tomultiplier (26.5 mm); this is explained by higher 
than average contamination in this device (but still 
within expected batch variability) . The vertical dis- 
tribution is also reproduced closely, except for the 
spike observed near 37 mm, the depth of the cath- 
ode wire grid. This is due to a previous contami- 
nation with radon arising in the xenon purification 
getters used prior to the first run. We subsequently 
found these to emanate ~1 atom/s, leading to sig- 
nificant implantation of 210 Pb in some detector sur- 
faces. Notwithstanding, both energy and spatial 
distributions of electron recoils are well reproduced 
by the simulation within the fiducial volume. 



2.2. Plating of radon progeny on the cathode 

Although other primordial decay chains also in- 
clude radon isotopes, 222 Rn ( 238 U chain) is the 
longest lived (Ti/ 2 =3.824 days). The issue of 
getter emanation had been previously identified 
in ZEPLIN-II [l(J, which relied on continuous 
gas recirculation through the same getter model 
(SAES PS11-MC500). Upon interruption of purifi- 
cation, the decay of the a activity was measured at 
Ti/2=3.83±0.1 days, supporting a 222 Rn contam- 
ination [l7j]. By then the same getters had been 
used to purify ZEPLIN-III xenon before the first 
run, which led to internal plating with 210 Pb. This 
did not pose a problem then, and we opted not to 
etch it away during the upgrade; we replaced the 
getter with a new model (SAES PS4-MT3). The 
effect of internal radon decay can be understood 
with reference to the chain sequence: 

222 Rn (a, 3.8d) -^ 21S Po (a, 3.1m) ^ 214 Pb (/3~ , 27m) ->• 
214 Bi (j8~ , 20m) -+ 214 Po (a, 164/tfs) -> 210 Pb (fi~ , 22y) -► 
210 Bi (fi~ , 5d) ^ 21 ° Po (a, 138d) -> 206 Pb (stable) (1) 

Once washed into the target the radon mixes uni- 
formly in the LXe. The first two a decays produce 
negatively charged atoms which promptly lose ex- 
cess electrons (the a-particle drifts to the cathode 
where it reduces to He). The next two (3~ decays 
create 214 Bi + and 214 Po + ions which cannot trans- 
fer their ionisation state to the rare gas and thus 
drift to the cathode. Finally, the a decay of 214 Po 
implants 210 Pb (T 1/2 =22.3 yr) within the wire grid, 
which bottlenecks the chain. This 210 Pb activity is 
still present on the cathode in the second run and 
has remained approximately stable at 19.4±0.3 de- 
cays/day over several months (per each of 210 Pb, 
210 Bi and 210 Po, as measured by the rate of 210 Po 
nuclear recoils). 

The electron-recoil background arising in the 
fiducial volume from this contamination is insignifi- 
cant in the context of the overall radioactivity bud- 
get (and was not included in the simulation); the 
same is true of nuclear recoils from the (a, n) reac- 
tion from 210 Po on the stainless steel wires. In fact, 
these populations can be useful as stable sources 
of ionisation; they allow, for example, the ionisa- 
tion yield for the deepest events to be monitored 
as a function of time. The Bateman decay equa- 
tions for the 210 Pb sub-chain indicate that the ra- 
tios of 210 Pb, 210 Bi and 210 Po become constant af- 
ter ~2 years (which is approximately the period 
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Table 1: Component-level analysis of electron recoil background and nuclear recoils from neutron elastic scattering in 
ZEPLfN-flf. The sum of simulated contributions (excluding components without measured contamination, indicated 
in brackets) is labelled 'SSR total'. The additional rejection efficiency provided by the veto is indicated: 'ptag' for 
prompt tags, which remove 7-rays, and 'dtag' for delayed tags, which reject neutrons [f| (weighted averages are 
indicated under 'SSR total'; these agree with the measured tagging efficiencies). After veto tagging, the neutron 
event expectation in a f year long dataset with typical signal acceptance is some fO times lower (~0.4 events). 



Material 


mass, kg 


e-recoil, dru^ 


ptag 


n/year-f 


dtag 


Krypton-85 




0.007 


-0 






Ceramic fecdthroughs 


0.9 


0.08 


0.30 


1.35 


0.58 


Photomultiplicrs 


4.2 


0.40 


0.26 


0.74 


0.58 


Rock (halite) 




-0 


-0 


0.53 


0.58 


Polypropylene shield 


1,266 


0.25 


0.04 


0.10 


0.58 


Scintillator modules 


1,057 


0.09 


-1 


0.03 


-1 


Copperware 


-400 


(<0.1) 


0.10 


(<0.15) 


0.58 


Lead castle 


-60,000 


0.01 


0.54 


-0 


0.58 


Radon-222 


1 m 3 


0.03 


0.19 


-0 




Muon- induced 








-0.3 


-1 


SSR total 




0.86±0.05 


0.28 


3.05±0.5 


0.58 


SSR data 




0.75±0.05 


0.28 


n/a 




(FSR [6] 




14.5±0.5 


(36±18)* 


") 



' events/kg/day/keVee at 10 keVee 

* single scatters in 2,370 kg-days over 5-50 keVnr (unity detection efficiency) 

* FSR dataset expectation was 1.2±0.6 events in WIMP search box (net 128 kg-days) 



elapsed between the start of the two runs); when 
this so-called transient equilibrium is reached, the 
relative spectrum of decay products becomes con- 
stant, and so does the mean ionisation yield per 
event. Monitoring the S2/S1 ratio for these low- 
energy events allows, for example, the accuracy of 
the electron lifetime correction to be checked. 

2.3. Photomultiplier 7-roys 

The original phototubes (2-inch ETEL D730Q) 
were purchased a decade ago and had relatively 
high contamination levels (250 ppb U, 290 ppb Th 
and 1350 ppm K), releasing 1,400 mBq per unit 
in 7-ray activity. This dominated by far the back- 
ground budget in the first run (14.5 dru). For this 
reason a second run of the experiment with up- 
graded devices was planned from the outset. A 
low-background product (D766Q) was developed 
by ETEL in collaboration with the project to en- 
able pin-by-pin compatibility with the existing ar- 
ray. All materials and sub-components were radio- 
assayed with HPGe and/or ICP-MS. The best pho- 
totubes have 35 mBq in 7-rays, which was achieved 
through a complete redesign of the device. Most of 
this activity is concentrated at the rear of the en- 
velope, away from the fiducial volume. Simulations 



of background contribution to ZEPLIN-III take the 
spatial distribution of activity into consideration. 
A value of 0.40 dru is predicted from internal 7- 
rays, which represents a 40-fold improvement over 
the previous array and makes this component only 
marginally dominant. The equivalent simulations 
for the first run reproduced data well, which gives 
added confidence in this calculation. Backgrounds 
arising from the j3~ activity internal to the pho- 
totubes, which are often neglected in this type of 
study, are discussed in Section |3~T1 

2.4- Krypton-85 and other intrinsic backgrounds 

Several LXe-intrinsic backgrounds can affect this 
and similar experiments. 210 Pb from 238 U and 
222 Rn chains and 228 Ra from the 232 Th chain, for 
example, undergo low-energy (3~ decay occasionally 
lacking significant x- and 7-rays to help flag these 
events ('naked betas'). Tritium (T!/2=12.3 yr, 
/?max=19 keV) is produced cosmogcnically and, on 
some accounts, it could generate as much as —1 dru 
if the xenon is stored on the surface for long peri- 
ods [3]. However, it is reasonable to assume that 
these contaminants can be removed effectively dur- 
ing xenon purification and, so long as radon em- 
anation into the detector can be controlled, they 
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should not come to dominate the background at 
low energy. We see no evidence for a substantial 
presence of these in the second run of ZEPLIN-III. 

More significant is the threat posed by 85 Kr, 
long recognised as a serious background for WIMP 
searches. 85 Kr is an anthropogenic (3~ emitter 
(T 1 / 2 =10.76 yr, /3 mox =687 keV); its presence in 
the atmosphere is mainly due to the reprocessing 
of spent nuclear fuel. The pre- nuclear 85 Kr/Kr ra- 
tio was as low as 3xl0 -18 in the early 1950s, but 
the present-day ratio is -~lxlO -11 fljj ]. 

The xenon used in ZEPLIN-III was sourced in 
the 1970s from an underground origin. An accu- 
rate measurement of its 85 Kr content was obtained 
in 1997 from a 5.9 g sample jicj], from which we de- 
rive an age-corrected 85 Kr/Kr ratio of 1.5xl0~ 12 . 
The overall Kr content was subsequently reduced 
by cryogenic distillation to ~50 ppb Kr (by weight), 
corresponding to 0.2 dru at low energy in the xenon 
target. This was, however, a conservative estimate 
and the present data suggest much lower contami- 
nation, which can be confirmed independently. 

A minor (3~ decay (0.434%, /3 mQ x=173 keV) is ac- 
companied by a 514 keV 7-ray from the 85 Rb 9/2 + 
level (T 1/2 =1.015 /xs) (U. This allows a f3 - 7 de- 
layed coincidence analysis, with the 7-ray searched 
in the veto scintillator or in the LXe target. A 6- 
month dataset revealed ^9 such coincidences in the 
target, conforming to the expected time delay dis- 
tribution. This translates to 4xl0~ 21 85 Kr/Xe ratio 
(in present-day xenon this would represent a very 
low 150 ppt Kr). The decay rate is 0.007±0.002 dru, 
which is insignificant in ZEPLIN-III. 

2.5. Ceramics and copperware 

Along with the photomultiplicrs, two materials 
could make up a significant fraction of the radioac- 
tivity budget of the experiment from the outset: the 
^400 kg of copper (from which all structural ele- 
ments, vessels, fastening hardware parts, etc., were 
manufactured) and ceramic feedthroughs. Other 
items were screened and their locations chosen to 
minimise their impact: stainless steel parts (e.g. 
vacuum ports), indium seals, LN2 in the 30-litre 
reservoir. 

For copper (OFHC C103), conservative upper 
limits of 0.5 ppb for U/Th contamination were ini- 
tially adopted, although lower values existed in our 
database for similar material [l5l ]. Simulation of 
this contribution indicated an electron recoil back- 
ground of 0.35 dru at low energy. However, the 
measured energy and spatial distributions show no 



evidence for such a significant rate, and we revised 
the upper limit down to 0.1 dru for the primordial 
radioactivity. 

Cosmogcnic activation also produces several 
radio-isotopes, of which 60 Co (Ti/2=5.3 yr) will 
retain the highest activity a few years after ex- 
posure to atmospheric muons ceases. We adopt 
a production rate for 60 Co in natural copper of 
~50 atoms/kg/day [l8j]. Since the exact exposure 
time on the surface is uncertain we assume secular 
equilibrium, with a subsequent cool-off period of 
4 years (preceding the second run after deployment 
underground) . Simulation of the resulting 7-ray ac- 
tivity within the copperware produces a negligible 
0.01 dru in the fiducial volume. 

There are nearly 100 ceramic feedthroughs for 
signal and high voltage connections, weighing 900 g 
in total; these are distributed among the xenon 
chamber baseplate and the vacuum vessel flange. 
They are alumina based and their average contam- 
ination was measured at 105 ppb U, 270 ppb Th 
and 880 ppm K. The 7-ray background is not severe 
(0.08 dru), but they dominate the neutron back- 
ground in the detector as discussed later. 

2.6. Polypropylene shield and veto detector 

The ZEPLIN-III veto consists of a tight-fitting 
polypropylene shield enveloping the main instru- 
ment, itself surrounded by 52 plastic scintilla- 
tor modules forming barrel and roof sections. 
The gadolinium-loaded hydrocarbon moderates 
and captures internal neutrons very efficiently; the 
Gd capture 7-rays are detected by the scintillators. 
The neutron detection efficiency (delayed tagging in 
a 70 fis window) is 58%; internal 7-rays are detected 
with an average 28% efficiency (prompt tagging in 
±0.2 /is window). The prompt tag adds a further 
2% to the neutron efficiency. A detailed discussion 
of the hardware and physics performance is given 
elsewhere 0, [B[ ; there, we analyse the backgrounds 
contributed to the WIMP target (a prime design 
consideration) as well as background measurements 
in the veto detector itself. A tonne of scintillator 
within the lead castle provides confirmation of the 
local radiation environment and we find that the 
veto measurements agree with the background mea- 
sured within ZEPLIN-III itself. 

The readout (photomultipliers, electronics, ca- 
bling, etc.) are external to the hydrocarbon shield- 
ing to mitigate neutrons, which is achieved very suc- 
cessfully. These also have a negligible contribution 
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to the 7-ray background in the LXe target. How- 
ever, owing to their large mass, the polypropylene 
and the plastic scintillator do contribute apprecia- 
bly to the ZEPLIN-III budget. Contamination levels 
have been measured for the polypropylene (1 ppb 
U, 1 ppb Th and 5 ppm K) and the plastic scin- 
tillator (0.2 ppb U, 0.1 ppb Th and 0.2 ppm K). 
These values are close to the sensitivity of the HPGe 
and ICP-MS measurements of those samples, but 
this radiological contamination is entirely consis- 
tent with the interaction rate recorded in the veto 
itself. Furthermore, these have distinct energy and 
spatial distributions in ZEPLIN-III, as shown in 
Figure^ which confirm these contributions: the en- 
ergy spectra are generally softer and they increase 
the event rate especially at the top of the xenon tar- 
get, in contrast with other dominant components, 
such as the phototubes and ceramics, which are lo- 
cated below the xenon volume. 

We note that radioactivity in the plastic scintil- 
lator does not contribute to the overall background 
since the Q-value of the radioactive decay, most of 
which is deposited locally, is nearly always suffi- 
cient to trigger the veto (even if the highest energy 
7-rays are able to escape). On the other hand, the 
polypropylene generates 0.25 dru in the xenon tar- 
get and has a low veto probability of 4%. 

2.7. j-rays inside lead castle 

The 7-ray flux inside the empty ZEPLIN-III cas- 
tle includes external contributions from the labo- 
ratory walls, natural as well as cosmogenic activ- 
ity in the shield, air-borne radon and its metallic 
progeny plated out on the inner castle walls. 7- 
rays from Boulby rock (halite) are attenuated by a 
factor of '--TO 5 by the 20-cm thick lead shield [22J 
and are insignificant as a background in the WIMP 
target. The activity of the 222 Rn chain measured 
with a commercial monitor inside the empty castle 
was 2.4 Bq/m 3 , which is typical at Boulby. Simu- 
lation of 7-rays from this activity leads to 0.03 dru 
in ZEPLIN-III. Further in situ measurements with 
a HPGe detector place an upper limit on the back- 
ground from the castle (primordial and cosmogenic 
activity in the bulk shield and inner walls). The 
castle walls and roof are made from lead smelted 
into thin stainless steel containers to make 1.2- 
tonne chevron-shaped blocks; the lead had been 
underground for some two decades and is known 
to have low 210 Pb content (it was used to shield 
previous dark matter experiments at Boulby) . The 
shield flooring is OFHC copper, which was moved 



underground over a decade ago as 'new copper'. 
The HPGe spectra show only a vestigial amount 
of 60 Co. The combined contribution of the shield 
to ZEPLIN-III is -0.01 dru. 

3. Rarer types of electron recoil event 

3.1. (3- and ^-induced photomultiplier backgrounds 
A significant amount of potassium is involved 
in sensitising internal photomultiplier surfaces, 
namely the bialkali photocathodes. The 40 K thus 
located close to the fiducial volume can generate a 
response in the detector in several ways, with the 
1,461 keV EC 7-rays being perhaps the most ob- 
vious. However, the 89% branching ratio j3~ de- 
cay (/3 m(M ;=l,311 keV) can be as problematic (this 
represents 30 mBq in the D766Q). VUV-sensitive 
phototubes have thin quartz windows (—1 mm), al- 
lowing some internal /3-particles to reach the LXe. 
Others will generate bremsstrahlung photons in the 
quartz. Cherenkov photons are also produced in the 
quartz and can be detected in coincidence with the 
prompt scintillation. Other luminescence processes 
arise internally in the phototubes, such as 'dynodc 
glow' from electron impact, which is signal coinci- 
dent. Besides contributing some 10-20% to the de- 
tector trigger rate, these /3- induced events can lead 
to backgrounds which are harder to discriminate 
than single scatters in the fiducial volume. These 
might combine, for example, bremsstrahlung pho- 
tons interacting in the main volume with Cherenkov 
light from the quartz, thus decreasing the observed 
S2/S1 ratio of the fiducial interaction. Cherenkov 
emission from Compton electrons created by 7-ray 
backgrounds may be even more damaging as fewer 
photoclectrons may be involved in this instance. 

Figure [3] illustrates some of these processes; it 
shows Monte Carlo data for 40 K decays generated 
internally in the central phototube in ZEPLIN-III. 
The quartz windows are curved and <1 mm thick. 
The blue histogram (cumulative) represents energy 
deposits everywhere in LXe whereas the red his- 
togram refers to the active region above the cathode 
only (mainly 7-rays). Transmitted /3-particles and 
bremsstrahlung in the quartz window clearly dom- 
inate at low energies for interactions just above the 
phototube. In the fiducial volume, some 10 mm 
away, these contributions become negligible. 

As mentioned above, Cherenkov light generated 
inside the quartz may also contribute to SI, yield- 
ing potentially tens of photoelectrons. For exam- 
ple, a 1 MeV electron - e.g. a 40 K /3-particle or a 
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Figure 3: Energy spectrum in LXe from simulated 40 K de- 
cays internal to central photomultiplier; curves are scaled 
(x31) for correct fiducial rate from this contribution. The 
solid histograms are cumulative and represent interactions 
above the phototube window; these include bremsstrahlung 
photons generated in the quartz, transmitted /3-particles and 
non-brcmsstrahlung 7-rays. The spectrum above the cath- 
ode is shown in red, mostly due to 1,461 keV 7-rays. 



Compton electron from the 1,461 McV 7-ray - pro- 
duces ^100 photons/mm of quartz within the spec- 
tral response of these photomultipliers; some gen- 
erate photoelectrons even without leaving the win- 
dow, while others suffer total internal reflection at 
the liquid surface and strike the array again; there- 
fore, the photoelectron yield can be significantly 
higher than suggested by the nominal QE. This op- 
tical signal can both increase the size of a time- 
coincident scintillation pulse as well as reduce its 
apparent time constant, threatening single-phase 
(scintillation only) as well as two-phase detectors. 
We see evidence of Cherenkov emission in ZEPLIN- 
III. This is illustrated in Figure 21 where the pulse 
timing parameter r, which measures the mean ar- 
rival time of the SI signal, is shown for hducial scin- 
tillation pulses and for sub-cathode events (no S2) 
with the same range of apparent energies (30-60 
keVee). The latter are mainly due to /3-particles 
transmitted through the phototube windows - but 
a faster population is also clearly seen which we 
attribute to Cherenkov emission. We can rule out 
an explanation based on nuclear recoils on the un- 
derside of the cathode grid, which would indeed be 
faster than electron recoils and yield no ionisation, 
by comparing with the corresponding timing spec- 
trum for the upper grid surface. We conclude that 



pulses in this population are indeed faster than nu- 
clear recoil scintillation. More extreme examples 
can be identified if signals triggering in a single 
channel are included (such as the blue pulse shown 
in inset), but in this instance it is harder to prove 
an optical, rather than electrical, origin. 

Such Cherenkov processes can be dangerous if as- 
sociated with low-energy electron recoils in the fidu- 
cial volume. Besides Compton electrons generated 
by background 7-rays, this could involve, more gen- 
erally, j3~ decays with coincident 7-rays (many ex- 
amples exist in the U and Th chains). In ZEPLIN- 
III, the timing distribution of fiducial events with 
higher than average scintillation (lower S2/S1 ra- 
tio) at energies above 30 keVee does not correlate 
with shorter than average pulses; rather, it is fully 
compatible with that of electron recoils - although 
this conclusion stems from a relatively small expo- 
sure. In any case, this effect should be carefully 
assessed in future noble liquid instruments where 
scintillation is detected by photomultipliers. 

For this and other reasons optical cross-talk be- 
tween phototubes must be avoided, so that these 
effects remain localised to a single channel and are 
therefore easier to identify. It is also important to 
reduce the amount of LXe allowed around their en- 
velopes to minimise scintillation from internal /3- 
particles, and to prevent those photons from reach- 
ing the photocathode either by external or inter- 
nal paths. The copper structures which house the 
ZEPLIN-III phototubes allow a thicker LXe layer 
around their envelopes than with the original ones, 
since the upgrade photomultipliers arc slightly nar- 
rower. However, there is no significant direct path 
into the main volume and the quartz envelopes were 
metallised to prevent light from reaching the pho- 
tocathode directly. 

3.2. Multiple j-ray interactions: MSSI events 

Multiple interactions where one vertex occurs 
in a detector region which is optically coupled to 
the phototube array but from which no charge 
can be extracted can lead to a challenging back- 
ground. Multiple-scintillation single-ionisation 
(MSSI) events produce lower S2/S1 ratios than sin- 
gle fiducial scatters since the multiple SI pulses are 
time-coincident but there is only one S2. This was 
identified early on as a background mechanism in 
two-phase detectors. In ZEPLIN-III an initial as- 
sessment during construction showed that efficient 
rejection based on SI light pattern and S1-S2 ver- 
tex consistency would be necessary Q . In the first 
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Figure 4: Distribution of photoelectron mean arrival times 
(r) in Sl-like pulses for fiducial interactions (yellow his- 
togram) and for Sl-only events occurring between the pho- 
tomultiplicr windows and the cathode grid. The apparent 
energies are 30—60 keVee for both histograms. A popula- 
tion of fast pulses is clearly visible, which we attribute to 
Cherenkov photons generated in the phototube windows by 
internal /3 — emission. Inset: array time response for three 
events with 30 keVee reconstructed energy (approximately 
55 photoelcctrons); the red trace is a fiducial event with 
r=35 ns; the black trace, probably a scintillation event con- 
taminated by significant Cherenkov light, is chosen from the 
black histogram with r=13 ns. The pulse in blue has t=7 ns 
and occurs in only one channel (the other cases both trigger 
7 channels). 



run this event topology became perhaps the most 
challenging, in particular because these populations 
cannot be reproduced accurately with calibration 
point sources. 

A detailed calculation is difficult when the light 
collection for the dead vertex is not known pre- 
cisely. In ZEPLIN-III, this is the case in two re- 
gions which yield no charge: i) the reverse field re- 
gion below the cathode grid, where light collection 
is high but variable due to the proximity to the 
photomultiplicr windows; ii) the peripheral region 
near the side walls, where the light collection is low 
and also uncertain. In any case, it is instructive 
to assess what effect these events will have on dis- 
crimination if not successfully rejected. Consider 
two 7-ray interactions, one in the fiducial region 
and one in a dead region, which we assume to have 
equal light collection and a constant spectrum of en- 



ergy deposited (as suggested by the Klein-Nishina 
formula for low-angle scattering of high energy 7- 
rays). The energy is reconstructed from the coin- 
cident SI components, but there is only one S2. 
The number of such double scatters adding up to 
energy E is P(E)dE = P^EdE, where P is the 
single-scatter probability per unit energy. So, this 
effect becomes more frequent with increasing recon- 
structed energy, as might be expected. The discrim- 
ination parameter, S2/S1, will suffer most when the 
'good' vertex produces the smallest detectable S2 
signal (E t h) and SI takes its apparent energy from 
the dead vertex only (S2/Sl>E t h/E). However, 
the (x,y) positions recovered from SI and S2 in- 
dependently are likely to be inconsistent for these 
events, leading to their probable rejection. There- 
fore, the good SI vertex must produce a minimum 
viable signature and, in practice, E t h ~ 2.5 keVee, 
the SI energy threshold. This trend is indeed ob- 
served in our data before these events are removed 
by the vertex reconstruction algorithm. A sim- 
ilar calculation can be done for neutron interac- 
tions. Assuming, very approximately, an expo- 
nential recoil spectrum with characteristic energy 
Eq, then the spectrum of double- vertex events is 
P(E)dE(xEexp(—E/E ) dE (i.e. the relative num- 
ber of such events still increases with E). 

4. Radioactivity neutrons 

The agreement between data and simulations 
obtained for the spectrum of electron recoils at 
low energy confirms the radiological contamina- 
tions adopted for each component; these were sub- 
sequently used to derive neutron production rates in 
those materials. Muon-induced showers, in partic- 
ular photo-production in electromagnetic cascades, 
are the only noteworthy source of non-radioactivity 
neutrons. This contribution is not dominant in a 
detector of the size of ZEPLIN-III (H, HU and, in 
addition, these neutrons are vetoed effectively. The 
dominant sources of fast neutrons are, therefore, 
spontaneous fission (mostly of 238 U) and (a,n) re- 
actions arising from the U and Th chainsjj 

Absolute neutron spectra were calculated with 
SOURCES-4A and -AC 
sections extended to 10 MeV 



26( with (a,n) cross- 
as described in 



2 An a activity of 10 Bq is present in the polypropylene 
due to 152 Gd, but this is not a significant source of neutrons 
(£ Q =2.8 MeV). 
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Ref. [13 and complemented with EMPIRE2.19 [H| 
cross-sections and branching ratios for transitions 
to excited states [29[. The (a,n) reaction domi- 
nates in light materials, producing, e.g., nearly all 
neutrons emitted by the Boulby rock (considered 
pure NaCl here), very few in copper (<10%) and 
practically none in lead. Yields in pure elemental 
materials are generally in good agreement with ex- 
perimental data, but it should be noted that the 
exact composition, uniformity and granularity of 
mixtures in detector components is a cause of un- 
certainty due to the very short range of a-particles. 

Primary neutron spectra are generated isotropi- 
cally from an appropriate position in the GEANT4 
simulation geometry and tracked to the LXe tar- 
get. In the case of the photomultipliers this is done 
at sub-component level, since some materials are 
further away from the fiducial volume than oth- 
ers. Each photomultiplier emits just over 2 neu- 
trons per year; this is dominated by a small amount 
of borosilicate glass still present in this model; al- 
though there is scope for further improvement, this 
value represents a 50-fold reduction relative to those 
used in the first phase of the experiment. 

The most challenging calculation is the contri- 
bution from the laboratory rock. It is computa- 
tionally intensive and requires accurate modelling 
of the shield: in spite of the high overall attenuation 
factor of ~10 5 , any gaps for pipework, cabling and 
alignment tolerances need to be taken into account. 
To ensure an accurate result a detailed CAD solid 
model of the veto and shield was imported into the 
GEANT4 model. The Boulby rock has 65 ppb U 
and 130 ppb Th assumed in secular equilibrium [30| . 
Neutrons are generated in the rock to a depth of 
3 m and tallied initially on a test surface in the lab- 
oratory. The simulated integral flux is adjusted at 
this stage to measurements at Boulby [3l[, includ- 
ing re-entrant (i.e. backscattered) neutrons. The 
differential spectrum is released with a cosine bias 
from a generator sphere containing the experiment 
to generate an isotropic and uniform flux within it. 

Yearly rates of single scatters in 5-50 keVnr in 
the fiducial volume and the component-level break- 
down are given in Table I and shown in Figure [SJ 
The total of 3.05 events per year is dominated by 
the ceramic feedthroughs (44%), followed by the 
photomultipliers and rock neutrons. We note that 
the background expectation predicted for a one- 
year dataset with realistic signal acceptance is con- 
siderably lower; typical signal detection efficiencies 
(including 90% operational duty cycle, 50% accep- 



tance in S2/S1 discrimination parameter, a further 
25% loss from software selection cuts), bring this 
down to 1.0 n/year; the rate in anti-coincidcncc 
with the veto is 60% lower, at 0.4 n/year. Finally, 
the detection threshold is likely to be a little higher 
than 5 keVnr, so a value just below 0.3 n/year is ex- 
pected. This is very close to the original goal of the 
experiment. The veto-coincident data analysed so 
far contains zero delayed coincidences (dtag) within 
the signal acceptance region, but this constrains the 
latter event rate only to <3 n/year at 90% CL (note 
that vetoed events have not been subject to a blind 
analysis since they contains no information on the 
presence of a signal). 




nuclear recoil energy, keVnr 



Figure 5: Simulated differential and integral rates of neutron 
single elastic scattering in the fiducial volume of ZEPLIN-III. 

The event rates discussed above are for single 
scatters in the 6.5 kg fiducial volume with no ad- 
ditional interactions (elastic or otherwise) in the 
12.5 kg active volume. Where more than one elastic 
scatter occurs, the low S2/S1 ratios (and short scin- 
tillation time constant) characteristic of nuclear re- 
coils are essentially preserved. Although the event 
can still be rejected from the multiple S2 pulses, it 
is often useful to establish the distribution of scat- 
tering multiplicity. This may provide additional 
constraints to the expected rate of single scatters, 
which arc an irreducible background for WIMP 
searches. In ZEPLIN-III, some 50% of events with 
energies reconstructed in the 5-50 keVnr region 
are single interactions. The mean scattering length 
varies quite significantly with neutron energy, but 
it is much less sensitive as a function of recoil en- 
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ergy (A ~ 30 cm). For systems much larger than 
this value, self-shielding can render the fraction of 
single scatters very small. Multiple scatters can be 
identified by additional S2 pulses in two-phase sys- 
tems (although they will become MSSI events if the 
second vertex yields no charge). However, the time- 
of-flight for radioactivity neutrons should not be ne- 
glected as a valuable discriminant. The mean time 
between scatters (above threshold) is only ~15 ns 
in ZEPLIN-III, but this effect can be readily seen in 
the structure of SI pulses during neutron calibra- 
tion, which points to its usefulness in large single- 
phase detectors. 

5. Discussion 

In this study we analysed electron and nuclear 
recoil backgrounds of radiological origin in the 
ZEPLIN-III experiment now operating in its sec- 
ond phase at Boulby. Energy spectra and spa- 
tial distributions are presented for electron recoils 
in the fiducial volume. The low energy rate is 
0.75±0.05 evt/kg/day/keVee, which represents a 
20-fold improvement relative to pre-upgrade levels. 
A further 28% reduction is possible due to the new 
veto detector installed around the WIMP target. 
Significantly, these vetoed events can be studied 
without compromising a blind analysis of the signal 
region. 

The electronic background achieved in phase II 
meets the original goal of the ZEPLIN-III pro- 
gramme. It implies a leakage of 2-3 events/year 
for a discrimination power of 1:7,800 (the value 
achieved in the first science run). Whether this 
can be matched in the second run depends on the 
optical performance of the photomultipliers, which 
unfortunately is not ideal. 

A long campaign of component-level radio- 
assays (and Monte Carlo simulations to determine 
their contribution to the radioactivity budget) led 
to good matches to the measured distributions. 
We confirmed that the photomultipliers and the 
polypropylene shield dominate, while copper and 
85 Kr, identified early on as possible concerns, con- 
tribute very little. Similar studies for other experi- 
ments, namely XENON100 [HI, have been equally 
successful, confirming that this methodology pro- 
duces accurate results. 

The good agreement achieved for electron re- 
coils also confirms the primordial U/Th contamina- 
tion of key components, lending more certainty to 



the calculation of neutron backgrounds from nat- 
ural radioactivity. The total rate of single elastic 
scatters in the fiducial volume is 3 events/year in 
the energy range 5-50 keVnr, which translates to 
0.4 events/year in anti-coincidence with the veto 
for a realistic signal acceptance. This will not limit 
the sensitivity of the experiment. We also pointed 
out that neutron time-of-flight can be important to 
establish scatter multiplicity in large single-phase 
experiments. 

We examined event topologies which chal- 
lenge the high discrimination power which can 
be achieved for single-vertex electron recoils. 
Chcrcnkov light produced in photomultiplicr win- 
dows by internal j3~ activity or by Compton elec- 
trons from background 7-rays can be detected in 
coincidence with 7-ray fiducial scatters, thereby 
making scintillation pulses appear larger as well as 
faster - i.e. more similar to nuclear recoils. Dynodc 
glow is another cause of signal-coincident light, al- 
though it should affect both SI and S2 pulses. In 
this instance photon emission is delayed with re- 
spect to the original optical stimulus depending on 
the input optics configuration and its voltage bias. 
These mechanisms were identified several decades 
ago (see, e.g., [33| and references therein). In the 
context of dark matter searches, tests conducted at 
Imperial College in the early 1990s confirmed that 
dynode glow and Cherenkov light were observed 
when two photomultipliers face each other without 
an intervening scintillator [34j]; subsequently, the 
NAIAD experiment at Boulby documented a pop- 
ulation of fast noise events (r < 100 ns) which were 
attributed to dynode glow (see Fig. 1 in Ref . [35| ) . 
More recently, a WIMP search with CaF2 (Eu) scin- 
tillator at Kamioka Observatory found that the sen- 
sitivity below 10 keVee was limited by Cherenkov 
photons produced in light guides by Compton elec- 
trons caused by background 7-rays [36[ . This type 
of event may be difficult to identify if the light pat- 
tern is not too peaked in a single channel. We note 
that these effects should also concern single-phase 
experiments such as XMASS [37|, MiniCLEAN [H[ 
and DEAP-3600 [H]. 

We discussed also 7-ray MSSI events and how 
these degrade the discrimination power of two- 
phase detectors. They pose a real challenge to 
current experiments and their mitigation should 
attract significant design effort in future systems. 
Whereas single- vertex electron recoils have a steep 
distribution in the S2/S1 discrimination parame- 
ter, and their leaking into the nuclear recoil re- 
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gion is progressive, this is not the case for MSSI 
events. Fortunately, both the frequency and sever- 
ity of these events decreases for lower energies, 
where more signal is expected. 

Under the assumption that the discrimination 
power eventually observed is not hindered signifi- 
cantly by poor photomultiplier performance or in- 
complete rejection of MSSI events, ZEPLIN-III 
should achieve a sensitivity of ^lxl0~ 8 pb-year to 
the scalar WIMP-nucleon clastic cross-section. 
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